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ABSTRACT 

We conduct the first self-consistent numerical simulations of a recently discovered 
population of 47 large, faint (ultra-diffuse) galaxies, speculated to lie in the Coma 
cluster. With structural properties consistent with very large low surface brightness 
systems (i.e. /i(g,0)< 24 mag arcsec“^; re comparable to the Galaxy), the red colour 
((g-r) 0.8) and assumed low metallicity of these objects compels us to consider a sce¬ 
nario in which these are underdeveloped galaxies whose early {zc^2) accretion to an 
overdense environment quenched further growth. Our simulations demonstrate the effi¬ 
cacy of this scenario, with respect to available observational constraints, using progen¬ 
itor galaxy models derived from scaling relations, and idealised tidal/hydro dynamical 
models of the Coma cluster. The apparent ubiquity of these objects in Coma implies 
they constitute an important galaxy population; we accordingly discuss their prop¬ 
erties with respect to a ACDM cosmology, classical LSBs, and the role of baryonic 
physics in their early formation. 
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1 INTRODUCTION 

Low surface brightness galaxies (LSBs) remain an important 
test for the favoured ACDM cosmology, given the dominance 
of dark matter in these systems whose number density ex¬ 
ceeds that of normal galaxies (Bothun, Impey & McGaugh 
1997; Dalcanton et al. 1997). Although distinguished by 
classification from high-SB systems and exhibiting qual¬ 
itatively different halo properties at low radii (e.g. 
de Blok, McGaugh & Rubin 2001), their conformance to the 
Tully-Fisher relation supports a continuity with late-type 
discs (Zwaan et al. 1995; Schombert & McGaugh 2014). 
In the standard galaxy formation framework, these dif¬ 
fuse systems nominally formed within low initial density 
fluctuations (Fall & Efsthathiou 1980; Mo, Mao & White 
1998), resulting in blue, gas-rich and slowly evolving discs 
(McGaugh & Bothun 1994). The vulnerability of such ten¬ 
uous systems to external influence typically precludes their 
existence in high-density environments (Dekel & Silk 1986; 
Rosenbaum et al. 2009; Galaz et al. 2011). 

Improved imaging techniques reveal however a surpris¬ 
ing new population of such galaxies within rich cluster en¬ 
vironments. Specifically, we refer to the recent detection 
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in the Dragonfly Telephoto array of large diffuse galax¬ 
ies (van Dokkum et al. 2015a, VD15) in the Goma cluster, 
whose association with the cluster is proposed on the ba¬ 
sis of both their spatial distribution and a lower limit on 
their distance as implied by the unresolved nature of their 
stellar component. With further examination via high qual¬ 
ity GFHT imaging, and assuming a distance equivalent to 
Goma, the authors assert 47 of these objects as Ultra-diffuse 
galaxies (UDGs) with effective radii (re) ranging from 1.5 
to 4.5 kpc, and central surface brightnesses of 24 to 26 
mag arcsec“^ (Table 1). The membership of the largest 
of these objects in Goma has since been confirmed using 
spectroscopy from the Keck I telescope (van Dokkum et al. 
2015b). 

While classical blue LSBs have been previously detected 
in the Virgo cluster (Impey, Bothun & Malin 1988), the 
UDGs lie on the faint end of the red sequence of the (more 
massive) Goma cluster (Gavazzi et al. 2013), a dichotomy 
possibly related to environmental dependencies. The recent 
accretion and quenching of gas-rich LSBs in Goma is consis¬ 
tent with 1) the slow formation of LSBs within cosmologi- 
cally underdense regions; 2) their recent accretion to large- 
scale structures/filaments (Rosenbaum et al. 2009), and 3) a 
recent build-up of the faint end/an increase in the dwarf-to- 
giant ratio in the cluster red sequence since z = 0.2 (Lu et al. 
2009; Gavazzi et al. 2013). Moreover, other faint dwarves in 
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Coma display generally radial/anisotropic orbits, suggestive 
of recent accretion from the field or as part of subgroups 
(Adami et al. 2009). 

If applying stellar population models, however, the me¬ 
dian UDG colour (g-r) = 0.8 zb 0.1 is consistent with an old 
stellar disc passively evolving over a timescale up to z = 2, 
if assuming the typically low metallicity of diffuse LSBs (i.e. 
oxygen abundance 12+log(0/H)^8; McGaugh & Bothun 
1994). The paucity of quenched field galaxies within this 
mass range (Geha et al. 2012) thus implies a scenario in 
which the UDGs were quenched upon infall at high redshift, 
leading VD15 to speculate that these are examples of failed 
L*-type galaxies. 

This early quenching is supported by evidence of an al¬ 
most complete red sequence in a z = 1.8 analogue of Coma 
(JKCS041; Andreon et al. 2014). The apparent absence of 
UDGs within 300 kpc of the cluster centre, presumably via 
their destruction here, is also not consistent with their very 
recent accretion to a kinematically hot host. More gener¬ 
ally, the apparent ubiquity of these UDGs in Coma conflicts 
with the theoretical rarity of large diffuse LSBs such as Ma- 
lin I (Hoffman, Silk & Wyse 1992), motivating an alternate 
scenario for their origin. 

In this first theoretical study of these objects, we adopt 
self-consistent numerical methods for hypothetical UDG 
progenitors to ascertain their evolution within high-density 
environments modelled on the Coma cluster. Guided by ob¬ 
servational constraints from VD15, we devise a theoretical 
template for the scenario in which these objects are nor¬ 
mal galaxies accreted as satellites and quenched soon after 
the standard epoch of disc formation (z = 2). Section 2 de¬ 
scribes our UDG model, while Sections 3 and 4 describe its 
hydrodynamical and tidal interactions within a rich cluster. 
Section 5 concludes this study with a discussion of these 
results. 


2 A HYPOTHETICAL MODEL FOR A UDG 

PROGENITOR 

Our study adopts a numerical method introduced in 
Bekki & Couch (2011) and Bekki (2014). The N-body ideal¬ 
isation of the UDG progenitor is constrained by the observed 
properties (Table 1), in particular the estimated median stel¬ 
lar mass of 6 X10^ M©, inferred from their median (g—r) and 
a tight mass-colour relation revealed in the GAMA survey. 

For an initial stellar (disc) mass (Md) of 10^ M©, 
abundance matching in cosmological simulations (e.g. 
Munshi et al. 2013; Behroozi, Wechsler & Conroy 2013) 
suggest a dark matter (DM) halo mass (Mh, within the 
virial radius rvir) up to ~10^Md, although rotation curves 
obtained from observed samples suggest these simulations 
produce too much substructure, attaining instead a factor 
10 lower prediction (Miller et al. 2014). We adopt an inter¬ 
mediate value for our UDG progenitor, and assume a mean 
density comparable to the Galaxy (assuming rvir,MW = 258 
kpc and Mh,MW = 1x10^^ M© Klypin, Zhao & Somerville 
2002), giving a rvir of 83 kpc. We assume an exponential 
disc morphology, with a density defined for the cylindrical 
radius r and height z: 

dr, 2 ) oc exp (-^)sech^(^). 


Table 1. Summary of a) properties of observed UDGs b) our 
reference model parameters 


(a) UDG Parameters 

Mean (Range) 

/r(g,0)/mag arcsec“^ 

25 (24-26) 

re /kpc 

3 (1.5-4.5) 

Sersic index, n 

1 (0.5-1.5) 

Stellar mass, Mg/MQ 

6x10^ (1x10^-3x10®) 

Gluster-centric radius/kpc 

- (>300) 

(b) UDG Model Parameters 

Value 

N°. DM halo particles 

1x10® 

N°. Stellar particles 

3x10® 

N°. Gas particles 

2x10® 

DM halo mass, Mh (Mq) 

3.2xl0i® 

NFW concentration, cnfw 

5.0 

DM halo Tvir (kpc) 

83 

Stellar (disc) mass, M^ (M©) 

10® 

Stellar scalelength, r^ (kpc) 

1.7 

Stellar scaleheight, (kpc) 

0.34 

Gas mass. Mg (Mq) 

5x10® 

Gas scalelength, rg (kpc) 

5.1 

Gas scaleheight, Zg (kpc) 

0.34 



Figure 1. (Top) Initial rotation curve (assuming V(r) = {GM(< 
r)/r)°-^) for the UDG progenitor model; (bottom) simulation 
BOUND in the Vrot/B-band Magnitude phase space (black filled 
circles); its evolution (with direction highlighted with the red ar¬ 
row) is consistent with an observed LSB sample (black crosses) 
with best-fit curve (dashed line), reproduced from Ghung et al. 
( 2002 ). 

The scalelength ig is derived a size-mass" scal¬ 
ing relation (Dutton & van den Bosch 2009; 
Ichikawa, Kajisawa & Akhlaghi 2012) whose shallow 
slope at this stellar mass range (a ~ 0.15, with respect 
to simple self-similar models based on virial relations 
where a ~ 0.33) is believed to reflect the role of feedback 
(Mo, Mao k White 1998; Dekel k Woo 2003; Shen et al. 
2003). 

Assuming rgff = 1.67rd, the sizes of the observed UDGs 
can be attained if we thereafter apply a linear scaling (Ag) to 
ra. If the UDGs are assumed as classical LSBs, Ag is a proxy 
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for variations in the initial spin parameter; our adopted 
As = 1.5 lies at ~lcr of typical spin distribution among sim¬ 
ulated halos (Maccio et al. 2007), and is consistent with the 
tendency for high spin halos to host LSB galaxies. 

This selection of free parameters (Md, As) yield a 
central surface brightness consistent with those observed 
(24</i(B,0)/(mag arcsec“^)<26). We further note that the 
resulting ratio le/^vir lies within 1 cr of a tight linear relation 
identified by Kravtsov (2013), consistent with the galaxy 
size being set by the halo’s initial specific angular momen¬ 
tum (Mo, Mao & White 1998). 

We adopt the NFW density profile for the halo, defined 
as: 


p(r) oc (r/rh) ^(1 + (r/rn)^) \ 

where ih = rvir/cNEW and cnfw is a concentration 
factor, here set to 5 in accordance with mass-redshift- 
dependent relations for the adopted Md (i.e. Maccio et al. 
2007; Munoz-Cuartas et al. 2011). This is similar to the 
characteristic cnfw among a sample of LSBs obtained 
from fits to high quality Ha-Hj (cnfw — 5 — 6; 
McGaugh, Barker & de Blok 2003), although the authors 
argue that a cored pseudo-isothermal density profile pro¬ 
vides a better fit than NFW. 

The interstellar medium (ISM) is considered isothermal 
(temperature 10^ K), and modelled with smoothed parti¬ 
cle hydrodynamics. We assume again an exponential density 
profile, as for the disc, and adopt a scalelength rg = 2.6rd 
(i.e. the sample mean from Kravtsov 2013), with a total mass 
(Mg) five times that of the disc, as inferred from scaling re¬ 
lations identified by Popping, Behroozi & Peeples (2015). 

Figure 1 illustrates the initial rotation curve of our 
model (with parameters summarised in Table 1). We do 
not readily find other curves of comparable extension and 
Vrot in the literature, noting also that, similar to previ¬ 
ous simulations of a similar nature (e.g. Kazantzidis et al. 
2011), Vrot will decline significantly due to substan¬ 
tial mass stripping. The curve is qualitatively similar to 
LSBs, in which the stellar disc is dynamically insignificant 
(Bothun, Impey & McGaugh 1997). Our model also lies on 
the same Tully-Fisher relation exhibited in the sample of 
LSBs compliled by Ghung et al. (2002). 

Following the method of Yozin V Bekki (2014), our 
SF/feedback model (see Bekki 2014, for more detail) is 
selected principally to avoid gas clumping in the gas-rich 
ISM while providing the low SF rates of generic LSBs 
(Schombert, McGaugh V Maciel 2013). To summarise, star 
formation occurs in the event of a convergent, cool region 
of gas with local density threshold >10°'^ cm“^, while the 
thermal component of Supernova feedback (0.9 x 10^^ ergs) is 
injected into the ISM over an adiabatic expansion timescale 
of 10® yrs. 


3 RAM PRESSURE STRIPPING IN A 
COMA-ANALOGUE CLUSTER 

In this section, we test the hypothesis that the UDG pro¬ 
genitor model can be efficiently quenched by ram pres¬ 
sure stripping (RPS) upon first infall and interaction 
with the intracluster medium (IGM) of a predecesor of 


Redshift, z 



Figure 2. Coma virial mass and concentration, cnfw? cls 
a function of lookback time/z, as derived from weak leas¬ 
ing methods (Gavazzi et al. 2009) and cosmological simulations 
(Munoz-Cuartas et al. 2011). 



X (Mpc) 


Figure 3. Equal-scale schematics of Coma-cluster models con¬ 
sidered in this study, at z = 0.023 and ^ = 2 (10 Gyr ago). 
White dotted lines represent characteristic NFW and virial radii. 
BOUND and infall orbits are conveyed with dashed lines, overlaid 
in the z = 0.023 case, on an optical image of Coma with green 
dots showing the locations of observed UDGs (reproduced from 
van Dokkum et al. 2015a). 


Coma aX z = 2. We derive a model of Coma using re¬ 
cent weak leasing measurements from deep CFHT images, 
which reveal a dynamical mass Mcoma (at z = 0.024) of 
5.1 xlO^"^ h“^Mo (Gavazzi et al. 2009). Assuming Ho = 70 
kms“^ Mpc“^, this estimate is low compared to previ¬ 
ous studies (i.e. Kubo et al. 2007) and is thus conserva¬ 
tive within the context of our work. We adopt statisti¬ 
cal measures of halo mergers/assembly, followed in cosmo¬ 
logical simulations (Fakhouri, Ma & Boylan-Kolchin 2010; 
Munoz-Cuartas et al. 2011; Ludlow et al. 2013) to attain 
the evolution of mass and cnfw as a function of redshift, 
with respect to their present values (Figure 2). 

To model the impact of hydrodynamic interactions on 
our model, we adopt the method of Bekki (2014), in which 
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a cluster satellite lies within a cubic lattice (with periodic 
boundary conditions) representing its local ICM. The SPH 
particles comprising this lattice have their density (and ve¬ 
locity, relative to the fixed satellite) vary according to a 
cluster-centric radial dependence described by a NFW pro¬ 
file (with cnfw as inferred above). For the extrapolated 
Mcoma at z = 2 , we obtain a total baryonic mass Micm 
from Lin et al. (2012), who use X-ray emission to establish 
0.1<MicM/Mcoma<0.15, and an ICM kinetic temperature 
interpolated from the cluster sample of Matsumoto & Tsuru 
( 2000 ). 

We assume the UDG progenitor is accreted onto the 
cluster with the cosmologically most common trajectory 
upon first infall (Benson 2005). Described in terms of the 
radial/tangential velocity at rvir, this trajectory does not 
vary significantly since z = 2. To avoid setting our satellite 
in thermal equilibrium with the ICM, the orbit commences 
at 1.5rvir, and is shown schematically in Figure 3. 

We find this scenario satisfies the hypothesis that rapid 
quenching of the UDG progenitor can occur at high z; their 
red colour can be facilitated by the drop in the SF rate upon 
first infall (relative to an isolated counterpart, as illustrated 
in Figure 4) and consistent with other numerical studies 
(e.g. Cen, Pop & Bahcall 2014). This occurs in spite of the 
exclusion, in this simulation, of cluster tides that would oth¬ 
erwise diminish the halo restoring force. While not all gas 
exceeds escape speed, being retained in the massive halo, its 
gas surface density is insufficient to permit further SF (gen¬ 
erally lying below the typical threshold density for SF of 
~10 M 0 pc“^). The RP stripped ISM shows a qualitatively 
similar (jellyfish) morphology to Virgo cluster satellites un¬ 
dergoing a similar process (e.g. Chung et al. 2007). 

The uncertainty introduced by the extrapolation 
of Coma’s properties at higher z should be high¬ 
lighted, given for example, the uncharacteristically high 
ICM temperature of Coma for its X-ray luminosity 
(Pimbblet, Penny V Davies 2014). Similarly, a recent anal¬ 
ysis of the GIMIC cosmological simulations for quenched 
galaxies residing in (I 0 ^^“^^ M©) host halos illustrates a 
highly variable RP efficiency with redshift (with primary 
cause as yet undetermined; Bahe V McCarthy 2015). The 
authors note, however, that quenched galaxies preferentially 
lie in overdense regions of the ICM (by a factor 10, with re¬ 
spect to the mean at a given cluster-centric radius; see also 
Tonnesen V Bryan 2008) in which case our model, devoid of 
such substructure, can be deemed conservative. 

We further emphasise that our model is of similar mass 
to the dwarf-type model of Bekki (2014), which was also 
demonstrated to be efficiently quenched in a Coma-type 
cluster environment during first infall. That his dwarf-type 
model was HSB, as contrasted with our high spin/diffuse 
model, would further indicate that the quenching efficiency 
demonstrated here is not sensitive to our assumptions re¬ 
garding the construction of the UDG progenitor model. 


4 SURVIVAL OF UDG PROGENITORS IN 
CLUSTER TIDAL ENVIRONMENTS 

Our UDG progenitor is hypothesised to reside in Coma 
over several crossing timescales (xcr —I Mpc/IOOO kms“^~I 
Gyr). Previous studies have highlighted how successive peri- 
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Figure 4. (Top panels) Face-on gas surface mass density distri¬ 
bution for the UDG progenitor during first infall into a Coma- 
analogue at 2 ; = 2. The simulation time is denoted at the top 
left, with effective and optical radii shown with red dotted lines; 
(bottom) star formation history of the aforementioned simulation 
(red, with cluster-centric radius shown with the blue dashed line), 
compared against an isolated example (grey). 

centre passages or harassment by other satellites can trans¬ 
form (Sersic index) n = I discs to ear her-type morpholo¬ 
gies (Moore et al. 1996; Kazantzidis et al. 2011). In this sec¬ 
tion, we compare the time-evolution of our models structural 
properties (illustrated in the diagnostic of Figure 5) with ob¬ 
served UDG constraints (Table I). 

The group’s tidal influence is modelled explicitly at each 
simulation time-step with the acceleration imposed upon 
the constituent stellar/halo particles of the UDG progeni¬ 
tor by the cluster halo’s gravitational potential. We adopt 
the fixed spherically-symmetric NFW profile of Section 3 for 
the halo mass distribution, and derive the gravitational po¬ 
tential at the location of an individual progenitor particle 
with its cluster-centric radius. 

The trajectory of the progenitor within the cluster po¬ 
tential is defined by a selection of orbits/redshifts. Orbit 
’infall’ refers again to the cosmologically most common in¬ 
fall trajectory at z = 2, adopted in Section 3 (Benson 2005). 
Orbit ’bound’ adopts the 2 : = 0.023 properties of Coma, 
with a ratio of apocentre to pericentre (rperi) set to the me¬ 
dian value determined from large-scale simulations ( 6 ; e.g. 
Ghigna et al. 1998), and rperi fixed at 300 kpc, similar to 
the apparent limit identified by VDI5. We accordingly test 
if UDG progenitors are destroyed below this limit with orbit 
’bound (<300 kpc)’, in which rperi is 150 kpc. We assume 
that the UDG model would occupy these bound orbits sub¬ 
sequent to its quenching during the initial parabolic orbit 
’infall’, as a result of dynamical friction experienced through 
interaction with a live cluster halo (not modelled here). 

Figure 5 conveys how, for our initial stellar disc 
mass Md = 10 ^ M©, the models following orbits ’bound’ 
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Figure 5. A comparison of parameters from simulations (left to 
right) ’bound’, ’infall’, and ’bound (<300 kpc)’, against observa¬ 
tional constraints (Table 1). From top to bottom: Sersic index n, 
stellar mass within g-band limit M^, effective radius rg, axial ra¬ 
tio b/a, 5 f-band central surface brightness /r(g,0), and the ratio of 
the mass (within rt) to the mass required to avoid tidal pruning. 


and ’infair, remain within observed ranges for UDGs. 
This is in spite of tidal mass loss (estimated from 
the stellar mass exceeding escape speed) and the pas¬ 
sive fading of the stellar disc. We estimate the latter 
from stellar evolution models (Bruzual & Chariot 2003; 
Portinari, Sommer-Larsen & Tantalo 2004), assuming an 
initial stellar metallicity of 0.03 Z©, comparable to faint 
Local Group satellites at z = 2 (Leaman et al. 2013) and 
consistent with passive fading since z = 1 — 2 to provide the 
observed ((g-r)) colour (VD15). 

These models remain disc-like and extended, consis¬ 
tent with VD15 who obtain a good fits in GALFIT for 
stacked images of the UDGs. These, together with 

deeper imaging and a multicomponent fit of the largest UDG 
(van Dokkum et al. 2015b), exhibit low n in the range 0.5 
to 1.5. Unlike in our similar study, in which 10^ M© dwarf 
satellites develop strong asymmetries and bulges within a 
group environment (Yozin & Bekki, in prep.), the disc sta¬ 


bilisation provided by a dominant halo, together with early 
removal of the ISM, leads to a persistently late-type mor¬ 
phology. These models also tightly match the Tully-Fisher 
relation, established from a sample of LSBs by Chung et al. 
(2002), for the duration of their evolution (Figure 1). 

The deviations from a unity n are caused principally by 
elongation of the disc and/or formation of tidal arms. These 
features are quantified with an axis ratio (b/a) estimated 
with the /i-weighted mean axial ratio among ellipses fit to 
isophotes (between the satellites’ central /i and /Xg = 28 mag 
arcsec“^ at 0.25 mag arcsec“^ increments) of the smoothed 
face-on /x distribution (using a gaussian kernel that matches 
the ~0.12 kpc FHWM of the CFHT imaging). The simple 
morphology of the models (lacking strong bisymmetries like 
bars due to a low disc mass) means this simple method is 
sufficiently robust for our purposes, as illustrated in Figure 
6. For orbit ’bound’, which we deem conservative in light 
of successive rperi at 300 kpc, b/a falls to the lower end 
permitted by observations (Figure 7); by contrast, ’infall’ lies 
at the upper end. This metric alone therefore suggests the 
UDGs occupy orbits in between these extremes, although 
we do not account for the effects of their inclination. 

For orbit ’bound (<300 kpc)’, we can demonstrate the 
effective destruction of the UDG progenitor, denoted by sub¬ 
stantial mass loss and disruption to the previously exponen¬ 
tial disc, clearly illustrated in mock GFHT images that do 
not compare well with the observed counterparts (Figure 6). 
The systems do not, however, resemble the ultra-compact 
dwarves often found concentrated at their respective clus¬ 
ter centres (Drinkwater et al. 2003). Although occupying a 
similar magnitude range to the UDGs, the small sizes of 
UGDs (<100 pc) are speculated to reflect the loss of up to 
98 percent of their original luminosity via tidal stripping 
(Bekki, Gouch & Drinkwater 2001). 

The bottom row of Figure 5 shows how this satellite 
disruption is related to an insufficient mass (mt) enclosed 
within its tidal radius rt when interacting with the cluster 
tidal field at pericentre. For an orbital radius r, tides ex¬ 
erted by the cluster mass M enclosed within r can yield 
substantial pruning of the satellite if mt/(3M(rt/r)^)^l 
(Binney & Tremaine 2008). For orbits ’bound’ and ’infall’, 
we find this mass ratio of the order of unity for several Ter, 
consistent with their apparent robustness in morphology. In¬ 
cidentally, we note that VD15 adopt this same formulation 
to predict, for the UDGs estimated stellar mass, a stellar- 
halo mass ratio (Md/Mh) to account for these objects sur¬ 
vival at 300 kpc; their value of ~0.04 within rt is consistent 
with the initial conditions of our progenitor models (Section 
2 )- 


5 DISCUSSION AND CONCLUSIONS 

Using constraints from the first detections of large faint 
galaxies (UDGs) in the Goma cluster (VD15), we have de¬ 
veloped a theoretical template for their evolutionary history. 
The exceptionally diffuse stellar component would imply a 
late infall as deemed common theoretically and observation- 
ally among LSBs (Rosenbaum et al. 2009), yet the degener¬ 
acy associated with their red colour permits a passive fad¬ 
ing since as early as z = 1 — 2 (depending on the assumed 
metallicity). In this context, we have examined the intrigu- 
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Figure 6. (Top) CFHT imaging of UDGs, illustrating the range 
in size (reproduced from van Dokkum et al. 2015a); (bottom pan¬ 
els) mock CFHT images (^f-band surface brightness) of (from bot¬ 
tom to top) ’bound’, ’infall’, and ’bound (<300 kpc)’, viewed face- 
on at simulation time (T =) 2.1 Gyr (first pericentre passage), 
3.9, and 5.6 Gyr (second passage ). Coloured labels provide the 
instantaneous 6/a,, /i,(g,0), n and rg. 


ing possibility that these UDGs are an example of satelite- 
quenching at high redshift. 

While there are no previous examples of large red LSBs 
in clusters, the UDGs possess several characteristics (n, 
/i(g,0), b/a) similar to Local Group dSphs (VD15). It is in¬ 
triguing therefore that in an inventory of local resolved stel¬ 
lar populations, the Galactic companions typically formed 
30-50 percent of their mass (but in some cases, only 10- 
20 percent) by z = 2 (Weisz et al. 2011). This downsiz¬ 
ing (lower mass galaxies showing more prolonged formation 
epochs) supports the role of reionisation in these objects’ 
early formation, where gas heating by the UV-background 
suppresses SF by both preventing gas collapse/accretion and 
heating the ISM, delaying cooling until z = 1 (Babul & Rees 
1992; Skillman et al. 2003). The recent finding of a central 
depression in the surface brightness profile (relative to an 
exponential disc) of the largest UDG by van Dokkum et al. 
(2015b), a feature common amongst dSphs, leads them to 
speculate further the role of early stellar feedback in su- 


pressing SF at high-z and contributing to gas loss due to 
ram pressure stripping (e.g. Stinson et al. 2013). 

It remains to be demonstrated if these mechanisms 
would be efficient for our UDG progenitor (with Vmax of 40- 
50 kms“^), but we note that the construction of our model 
(Section 2), in particular the choice of As, accommodates 
the scenario in which the UDGs are underdeveloped HSB 
galaxies. Since log(As)^a, our progenitor model at z = 2 
can be assumed to have nominally developed a factor ~10 
more massive stellar disc if not quenched so early. 

We presently find tentative support for an early accre¬ 
tion/quenching scenario with the deprojection (with an Abel 
integral identity) of the UDGs spatial distribution (Figure 
7), which reveals a number density (as a function of ra¬ 
dius) that can be reasonably fit with a NFW profile with 
cnfw of 3, consistent with the satelite population at large 
(Lin, Mohr & Stanford 2004). 

If ignoring the small sample size, this implies a relaxed 
distribution, in a cluster whose two-body relaxation time is 
of order a Hubble time. However, the more centrally located 
dwarf satellites of Goma tend to be older (in spite of a wide 
range of stellar ages depending on the chosen mass proxy; 
Smith et al. 2012). Figure 7 indicates the 47 known UDGs 
exhibit no clear trends among their properties as a func¬ 
tion of cluter-centric radius, although VD15 note that the 
intracluster light may be obscuring some UDGs. 

Our tidal diruption model is simple, insofar as the tidal 
harassment from other Goma satellites is not explicitly in¬ 
corporated (and in which case our mass loss predictions may 
be lower limits). We have assumed at this stage that the 
UDGs, as low mass systems in a kinematically hot cluster 
(cr^lOOO kms“^; Danese, de Zotti & di Tullio 1980), are not 
signihcantly disrupted even on long timescales, as demon¬ 
strated in previous studies (e.g. Tormen, Diaferio & Syer 
1998). In general, dwarf galaxies are weakly inhuenced by 
dynamical friction in clusters (i.e. drag force oc satellite 
mass), with the associated orbit mostly avoiding regions of 
strong harassment (Smith, Davies & Nelson 2010). Diffuse 
systems such as the UDGs may remain quite vulnerable 
(Moore et al. 1996), but we have tentatively demonstrated 
here that the halo masses associated with these systems can 
provide some measure of stability. 

We cannot state however if this stability is suffi¬ 
cient in the event that the UDGs were hrst accreted as 
part of a group, as exhibited by simulated dwarf satel¬ 
lites (De Lucia et al. 2012). Our models have assumed a 
static host at the two epochs concerned {z = 0.023 and 
z = 2), yet Goma-cluster analogues within cosmological sim¬ 
ulations accrete >5 group-mass (10^^ M©) halos within the 
last Gyr alone (Fakhouri, Ma & Boylan-Kolchin 2010). The 
dynamical state of major subgroups in Goma (e.g. NGG 
4839) as pre- or post-merger is uncertain (Burns et al. 1994; 
Sanders et al. 2013), but the (radial, anisotropic) kinematics 
of other faint dwarves tentatively suggest a signihcant recent 
infall from the held or as part of subgroups (Adami et al. 
2009). 

Besides pre-processing in the more tidahy disruptive en¬ 
vironment of a group, tidal interactions between galaxies ac¬ 
creted as a group/hlament are also quite feasible (Gnedin 
2003). This scenario is not compatible with our hypoth¬ 
esized high-z quenching and \ow-z late-type morphology 
of the UDGs, but is consistent with cosmological simula- 
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Figure 7. (Top left) Deprojected number density of the 47 de¬ 
tected UDGs (using an Abel integral identity), normalised by the 
central density of a best-fit NFW profile with cnfw of 3 (see text; 
red dashed line), as a function of cluster-centric radius; (middle 
left) 5 f-band Magnitude, (middle right) b/a, (middle left) /r(g,0), 
and (middle right) Tq as a function of projected cluster-centric 
radius for the 47 detected UDGs, with a linear-fit to each data 
set (red dashed line). 


tions which rarely find large intact LSBs in clusters (e.g. 
Avila-Reese et ah 2005). 

To conclude, it is entirely possible that the ob¬ 
served UDGs are classical LSBs, recently accreted, and 
even on the verge of destruction. The recent deep imag¬ 
ing of the largest UDG evidence does not support its 
tidal disruption, although this may not be representative 
of the population (van Dokkum et al. 2015b). We might 
also speculate if they are the cluster-bound gas-stripped 
remnants of another population, recently detected in the 
ALFALFA survey, constituting exceptionally gas-rich (for 
their stellar mass) dark galaxies, typically with high spin 
parameters (Huang et al. 2012). The kinematics of the 
UDGs would certainly aid in identifying their relationship 
with their hosts’ dynamical history e.g. subgroup mergers 
(Vijayaraghavan, Gallagher & Ricker 2015). With this in¬ 
formation, and the dynamical mass/metallicities obtainable 
in future long-exposure LSB-optimised observations, we look 
forward to refining our theoretical model. 
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